arXiv: 1503.07514v2 [hep-th] 28 Aug 2015 


Hydrodynamics beyond the Gradient Expansion: Resurgence and 

Resummation 

Michal P. Heller 1,2 ’[] and Michal Spaliriski 2,3 ’[j] 

1 Perimeter Institute for Theoretical Physics, 

Waterloo, Ontario N2L 2Y5, Canada 
2 National Centre for Nuclear Research, Hoza 69, 00-681 Warsaw, Poland 
3 Physics Department, University of Bialystok, Lipowa 4C 15-424 Bialystok, Poland 

Consistent formulations of relativistic viscous hydrodynamics involve short-lived 
modes, leading to asymptotic rather than convergent gradient expansions. In this 
Letter we consider the Miiller-Israel-Stewart theory applied to a longitudinally ex¬ 
panding quark-gluon plasma system and identify hydrodynamics as a universal at¬ 
tractor without invoking the gradient expansion. We give strong evidence for the 
existence of this attractor and then show that it can be recovered from the divergent 
gradient expansion by Borel summation. This requires careful accounting for the 
short-lived modes which leads to an intricate mathematical structure known from 


the theory of resurgence. 
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Introduction. - The past 15 years have witnessed the rising practical importance of relativistic 
viscous hydrodynamics. One reason for this is the success of hydrodynamic modeling of 
quark-gluon plasma (QGP) in heavy ion collision experiments at the RHIC and the LHC 
and the realization that QGP viscosity provides a crucial probe of QCD physics [lj. Another 
motivation is the relation between black holes and fluids, which originated in the 1970s as 
the rather mysterious black hole membrane paradigm [2j. With the advent of holography, 
this connection has been promoted to a precise correspondence [3] shedding light both on 
the physics of QGP and liquids in general, as well as gravity. 

The perfect fluid approximation is widely used in astrophysics and this theoretical de¬ 
scription of relativistic inviscid fluids is rather well established [3]. On the the other hand, 
relativistic viscous hydrodynamics is much less well understood. One of the recent insights is 
to regard hydrodynamics as a systematic gradient expansion, much in the spirit of low-energy 
effective held theory [5]. 

However, the requirement of causality leads to a framework which necessarily incorpo¬ 
rates very large momenta (and frequencies). In all known examples this is accompanied by 
the appearance of short-lived excitations - nonhydrodynamic modes. It has recently been 
shown, in the context of AdS/CFT correspondence, that their presence leads to the diver¬ 
gence of the hydrodynamic gradient series for strongly coupled Af = 4 super Yang-Mills 
(SYM) plasma [6j. In view of this it is not clear whether or how a naive gradient expansion 
defines the theory. This is in fact a fundamental conceptual question concerning relativistic 
hydrodynamics as such. 

In this Letter we propose a definite answer: since the nonhydrodynamic modes decay 
exponentially, the system relaxes to an attractor regardless of when an initial condition is 
set. In the following we consider a simple situation in which this can be made completely 
explicit: the Miiller-Israel-Stewart (MIS) theory [5l [?, (8J specialized to a longitudinally 
expanding conformal fluid. We show that the attractor can be determined by relaxation 
from solutions which take the form of a transseries. The higher orders of this transseries are 
encoded in the divergent hydrodynamic gradient expansion, in line with expectations based 
on resurgence ideas [9]. 
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Miiller-Israel-Stewart theory.- The Landau-Lifschitz formulation of relativistic viscous hy¬ 
drodynamics [10] asserts that the evolution equations for the hydrodynamic fields - tem¬ 
perature T and flow velocity — are the conservation equations of the energy-momentum 
tensor 

(T^) =f«V + V(£)( rf v + vPu v ) + (1) 

where the shear stress tensor Yl liv is given by 

IF" = -r/cffi". (2) 

From a modern perspective, one could contemplate including on the RHS of Eq. (j2| all 
possible terms graded by the number of derivatives of T and u A If this is done to a finite 
order, as in Eq. ([2]), the resulting theory will not have a well-posed initial value problem due 
to superluminal signal propagation mm- 

MIS theory resolves this problem by promoting the shear stress tensor to an independent 
dynamical field which satisfies a relaxation-type equation: 

(T n u a d a + l)n^ = -r l a» I ' + ... , (3) 

where Tn is a phenomenological parameter (the relaxation time) and the ellipsis denotes 
several additional terms whose explicit form can be found in [5]. Linearization of the re¬ 
sulting theory is causal as long as Tr n > r\/s. This approach has enjoyed great success in 
describing the evolution of QGP [TS]. It has also been obtained as the long-wavelength effec¬ 
tive description of strongly coupled Af = 4 SYM plasma in the framework of the AdS/CFT 
correspondence oi r>.. 

The iterative solution of Eq. (J3]) generates the gradient expansion of the shear stress 
tensor leading to the appearance of an infinite number of terms on the RHS of Eq. ([2]). 
Their coefficients (the transport coefficients of all orders) are expressed in terms of ij and rn 
(and 4 more parameters in the general conformal case [5]). 

In the language of high energy physics, MIS theory can be regarded as an UV completion 
of Landau-Lifschitz theory. This is in contrast to the standard way of viewing it as a 
phenomenological model providing an effective description of some microscopic system only 
at late times. By treating MIS theory as an UV completion we mean that we consider a 
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hypothetical physical system such that the MIS theory describes it also at early times. This 
could be the case for systems where the approach to equilibrium is governed by a purely 
damped quasinormal mode. 


The setup.- To overcome the complexity of the MIS equations, we focus on the case of 
Bjorken flow [IS], which, due to a very high degree of symmetry, reduces them to a set of 
ordinary differential equations. The symmetry in question, boost invariance, can be taken 
to mean that in proper time-rapidity coordinates r, y (related to Minkowski coordinates t, z 
by t — Tcoshy and z = rsinhy), the energy density, flow velocity and shear stress tensor 
depend only on the proper time r. The MIS equations then take the simple form 


re = --e + 0, 

• Ar) Ai0 2 4 r n 0 
W = x- 


( 4 ) 


3r 2t] 2 3 r 

where the dot denotes a proper time derivative and </> = —11^, the single independent com¬ 
ponent of the shear stress tensor. The term involving Ai comes from the elided terms in 
Eq. (|3]); for details see j5j|- 

In a conformal theory, e ~ T 4 and the transport coefficients satisfy 


a 


r n = 


t n 


T ’ 


Ai =C X ,^, 


T] = C V S , 


( 5 ) 


where s is the entropy density and CVn, C\ 1 , C v are dimensionless constants. In the case of 
J\f — 4 SYM theory their values are known from fluid-gravity duality [3]: 


CVn — 


2 - log(2) 
27r 


= 2n' 


C '=4V 


( 6 ) 


The hydrodynamic attractor. - From Eq. Q one can derive a single second order equation 

for the energy density or, equivalently, the temperature: 

/ . \ 2 


T T\ 

tC t u~ + 3rCVn + ( 


, -L -L v_/ T 


3 T 


+ r)T + 


4 a 4 c. 


9 r 


+ 


rn 


9 r 


+ -r = a 


To simplify the presentation, we have set C^ = 0 in this equation as well as in Ec 
(10), and (13) below. 


( 7 ) 

(§, 
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To proceed further it is crucial to rewrite Eq. (JtJ) in first order form. Introducing the 
dimensionless variables w and / (as in Ref. im 


w = tT , 

the MIS evolution equation (J7|) takes the form 

C T nwff+4C rn f 2 +(w 


w 

f = T~, 

W 


( 8 ) 


i6 c. 


rn 


/ + 


4 C„ 16 C. 


rn 


2 w 

T = 0 ' 


( 9 ) 


and 0 together 


are 
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where the prime denotes a derivative with respect to w. Equations 
equivalent to Eq. ([7]) as long as the function w(t) is invertible. 

At large times (which translate to large w) we expect universal hydrodynamic be¬ 
haviour im- In phenomenological analysis of heavy ion experiments, usually based on MIS 
theory, hydrodynamic codes are initialized typically at w ~ 0.5, which corresponds roughly 
to a time r = 0.5 fm after the collision, with the temperature T = 350 MeV at the centre of 
the fireball at RHIC (see e.g. ro- Equation (J9]) indeed possesses a unique stable solution 
which can be presented as a series in powers of 1 jw\ 


f{w) = 3 + 9w 


4 C„ 8 c„a 


+ 




27w 2 






( 10 ) 


This is, in fact, the hydrodynamic gradient expansion. 

It is easy to see that linear perturbations around this formal solution decay exponentially 
on a time scale set by rn: 


5f(w) ~ exp 


2 a 


rn 


W W c rn 


1+ OR 

w 


( 11 ) 


This is precisely the short-lived mode introduced by the MIS prescription. In the language of 
the gravity dual to = 4 SYM theory this would be an analog of a quasinormal mode [6^ fl9] 
whose frequency is purely imaginary. 

The presence of this exponentially decaying mode suggests that Eq. g possesses an at¬ 
tractor solution. We propose that this attractor constitutes the definition of hydrodynamic 
behaviour. As discussed below, the presence of this attractor can be inferred without refer¬ 
ence to the gradient expansion, which, as shown in the following section, is in fact divergent. 
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FIG. 1. The blue lines are numerical solutions of Eq. ©> the for various initial conditions; the 
thick, magenta line is the numerically determined attractor. The red, dashed and green, dotted 
lines represents first and second order order hydrodynamics. The upper plot was made with 
parameter values appropriate for J\T = 4 SYM theory, while the lower plot has both rj/s and CVn 
increased by a factor of 3. Note that in the latter case the hydrodynamic attractor is attained at 
larger values of w, as expected. 

The existence of the hydrodynamic attractor is supported by examining the behaviour of 
generic solutions of Eq. ([8]), with initial conditions set at various values of w. 

As seen in Fig. [lj a generic solution rapidly decays to the attractor. Furthermore, the 
attractor appears to persist even at very small values of w, where hydrodynamics of finite 
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order becomes ill defined. Perhaps unsurprisingly, truncating Eq. (10) at first or second 


order gives results distinctly different from the attractor at very small w. The magnitude of 
this difference depends on the values of the transport coefficients. 

Assuming Af = 4 SYM parameter values, we see that adopting just the viscous hydro¬ 
dynamics constitutive relations provides a remarkably good approximation of the attractor 
for a wide range of w. In particular, this holds with an error smaller than 10% for w > 0.5. 

Examining the behaviour of / for w ~ 0 one finds two solutions, one of which is stable 


'2\JC r n + 

f{w) = - o /TV - + 0 \ W )• 


sVc. 


( 12 ) 


rn 


By setting the initial value of / at w ~ 0 arbitrarily close to the value given by Eq. (12), 
the attractor can be determined numerically with the result shown in Fig. [TJ 

Another way of characterizing the attractor is to expand Eq. ([£]) in derivatives of / - this 
is an analog of the slow-roll expansion in theories of inflation (see e.g. [20]). This leads to 
a kind of gradient expansion, but this is not the usual hydrodynamic expansion, since the 
generated approximations to / are not polynomials in 1/w. By choosing the correct branch 
of the square root which appears at leading order one can ensure that expanding the fc-th 
approximation in powers of 1/w one finds consistency with hydro at order k + 1. At leading 
order one finds 

9 . /fi4<A7TY 4- <W 

(13) 


, 2 w y/GlC/An + 9u; 2 

/ M = o - ^r- + -- 

o ou T n 


246Vn 

Continuing this to second order gives an analytic representation of the attractor which 
matches the numerically computed curve even for w as small as 0.1. 

Finally, one can also construct the attractor in an expansion around w = 0 starting with 


the f(w ) given by Eq. (12). It turns out that the radius of convergence of this series is finite. 
All three expansion schemes are consistent with the numerically determined attractor. 


Hydrodynamic gradient expansion at high orders.- In what follows we focus on the hydro- 
dynamic expansion, the expansion in powers of 1/w. It is straightforward to generate the 
gradient expansion up to essentially arbitrarily high order (in practice, we chose to stop at 
200). The coefficients f n of the series solution 

OO 

nw) = j2^ w ~ n 

n=0 


( 14 ) 











show factorial behaviour at large n, as seen in Fig. [2j This is analogous to the results 
obtained in Ref. [Bj for the case of AT = 4 SYM theory. 

In view of the divergence of the hydrodynamic expansion, we turn to the Borel summation 
technique. The Borel transform of / is given by 

fn 


/*«) = E ^ 


(15) 


71=0 


and results in a series which has a finite radius of convergence. Note that in Eq. (15) large 


w corresponds to small £. To invert the Borel transform, it is necessary to know the analytic 


continuation of the series (15), which we denote by /#(£). The inverse Borel transform 


Ir{w) = / die tf B {i/w)=w die f B (i) 

Jc Jc 


(16) 


where C denotes a contour in the complex plane connecting 0 and oo, is interpreted as a 


resummation of the original divergent series (14). To carry out the integration, it is essential 
to know the analytic structure of /#(£). 

We perform the analytic continuation using diagonal Pade approximants [21], given by 


the ratio of two polynomials of order 100. This function has a dense sequence of poles on 
the real axis, starting at £o = 7.21187, which signals the presence of a cut originating at 
that point |22j. This can be corroborated by applying the ratio method I2U. which allows 
the estimation of the location and order of the leading branch-cut singularity by examining 


the series coefficients. Specifically, if the function approximated by Eq. (15) has the leading 
singularity of the form (£o — £) 7 then for large n 


fn 


fo 


i + ^+0(4: 


(17) 


,/n+i n+1 V n + 1 ' n * 

Applying this formula one fords £o = 7.21181 (which is consistent with the pattern of poles) 
and 7 = 1.1449. 

As we argue in the Supplemental Material, the analytic structure of fs must involve 
further singularities on the real axis precisely in the following form: 


/*(«) = VO + Ko - 0 7 V0 + (20 - 0 27 V0 + ■ ■ ■ 


(18) 


where the functions h*,(£) are analytic and the ellipsis denotes further singularities at integer 
multiples of £ 0 - 
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Such branch-cut singularities of fs lead to ambiguities in the inverse Borel transform. 
Indeed, such a series is not Borel summable in the usual sense. It is, however, known that 
even in such cases a resummation is possible (see, e.g., Ref. (23]), but requires a nontrivial 
choice of integration contour. The freedom in the choice of integration contour leads to 
complex ambiguities 

poo 

Sfn(w) = e M ^w (19) 

Jk& 

where p is an odd integer reflecting the choice of Riemann sheet. For large w this becomes 

Sf R (w) « e iw ^r(fc 7 + l)/i*(fc£ 0 ) (ur 7 e- w& )* . (20) 


This ambiguity is a feature of the hydrodynamic series and its presence is an indication 
of physics outside the gradient expansion. We saw in the previous section that there are 
nonanalytic, exponentially suppressed corrections to the hydrodynamic series following from 
the presence of the nonhydrodynamic MIS mode. These have precisely the correct structure 
to eliminate the k — 1 ambiguity in inverting the Borel transform. Indeed, comparing Eq. 
(20) with Eq. ( [11] ) we are led to identify £ 0 with 3/2C r n and —7 with ( C n — 2C'y )/C' r n- 
Evaluating these combinations with parameter values appropriate for J\f = 4 SYM theory 
[Eq. (jo])] gives agreement to five significant digits. Both Eq. (20) and Eq. 0 receive 
corrections in 1 /w and we expect them to match also. 

The nonlinear structure of Eq. (|9]) suggests the presence of an infinite series of exponen¬ 


tial corrections, which are matched by further branch cuts in Eq. (18). In the following 
section we calculate these corrections and give strong evidence that they conspire to yield 
an unambiguous, finite and real answer for /r, up to a real constant of integration. 


Resurgence - The results presented so far suggest that Eq. ([ 9 ]) should possess a solution in 
the form of a transseries [23]: 

OO OO 

fW = ^ C m £l{w) m Y, am,nW~ n , ( 21 ) 

771=0 77=0 

where = w -7 exp(— w£, 0 ) while c and a m ^ n are coefficients to be determined by the equation. 
By direct substitution one can check up to high order that all of the coefficients a m ^ n in 


Eq. (21) are fixed uniquely apart from a^o, which can be absorbed into the constant c. 
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FIG. 2. The large order behaviour of the hydrodynamic series. The slope is consistent with location 
of the singularity nearest to the origin as given by Eq. ©• 


For each value of m in Eq. (21), the series over n is expected to be divergent - we have 
checked this for m < 2. Applying the Pade-Borel techniques discussed earlier leads to 
complex resummation ambiguities for each of these series. To obtain a meaningful answer, 
it must be possible to choose the single complex constant c in such a way that the result 
does not depend on the choice of integration contours and that the imaginary parts cancel. 

The key observation is that the ambiguity at the leading order of the transseries is propor¬ 
tional to fl, so it can only be canceled by terms of order m — 1 or higher. This cancellation 
determines the constant c 

c = r-e^ p r (7 + 1 ) / 4 0) (£o) (22) 


up to an arbitrary real number r, which is the expected integration constant for the first 
order differential equation ([9]). 

Resummation - Having provided strong evidence for the existence of an unambiguous and 
physically sensible result encoded in the transseries, we now invert the Borel sums for m < 2. 
For these calculations we used extended precision arithmetic (keeping one thousand digits). 
Inverting the Borel transform at each order of the transseries requires performing the in¬ 


tegration in Eq. (16). The analytic continuation by Pade approximants works well in regions 
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of the complex plane away from branch-cut singularities, so we take all of the integration 
contours to be straight lines at arg(£) = 7 t/4 . The integrals computed in this way are com¬ 


plex. The findings of the previous section suggest that by taking the sum as in Eq. (21) one 
should be able to choose the imaginary part of the constant c so that the result is real for 


some range of w. This is indeed the case and gives a value for Im(c) consistent with Eq. (22) 
(with p — — 1). A combined measure of error is the imaginary part of the result of the 
resummation - it remains very small (below 0.01% relative to the real part) for w > 0.25. 

We compared the result of the resummation with the numerically computed attractor, 
which required fitting the integration constant r = 0.049 [see Eq. (22)]. As seen in Fig. |3j the 
generalized Borel sum of the gradient series indeed follows the attractor. Note that to match 



FIG. 3. The hydrodynamic attractor (magenta), compared with the resummation result (cyan, 
dot-dashed) and the gradient expansion of order 1 (red, dashed) and 2 (green, dotted). 

the attractor we need to choose the coefficient r % —Re(c), while naively one might expect 
that the attractor should correspond to omitting the exponential terms. This suggests that 
the resummation of the gradient series contains exponentially decaying terms which are 
not canceled by the exponential terms from higher orders in the transseries. Consequently, 
instead of thinking about nonhydrodynamic modes in terms of perturbations about the 
gradient expansion one should more properly think of them as perturbations around the 
hydrodynamic attractor. 
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Conclusions - Recent advances in applying the theory of resurgence to quantum theories |25j- 
[30] have motivated us to rethink the foundational aspects of relativistic hydrodynamics. 
The root of the problem is that causality precludes us from regarding hydrodynamics as a 
truncated gradient expansion, yet the series itself is divergent. We propose to view hydrody¬ 
namics as an attractor which governs the late time behaviour of systems in their approach 
to equilibrium. In the context of boost-invariant flow in MIS theory we have constructed 
such an attractor in several ways which are all consistent with each other. 

By identifying the structure of singularities of the analytic continuation of the Borel 
transform of the hydrodynamic series in terms of nonhydrodynamic degrees of freedom, 
we hope that this Letter will provide a useful road map for understanding the meaning of 
higher order gradients in the case of AT = 4 SYM theory [6]. From the point of view of 
AdS/CFT, the exploratory studies described here suggest that the geometry constructed in 
the gradient expansion of fluid-gravity duality should be viewed as the leading term of a 
transseries containing the effects of quasinormal modes. 

From a phenomenological perspective, MIS theory includes explicitly transport coeffi¬ 
cients for terms up to second order in gradients, but it generates a gradient expansion to 
all orders. The transport coefficients can only match real QCD plasma up to second order. 
One may then wonder about the effects of all of the higher order terms EU which cannot 
be matched by MIS theory. Our hirelings suggest that the attractor which governs its late 
time behaviour is not very sensitive to the higher order terms (even when the gradient series 
is resummed). This makes it less surprising to learn that MIS theory can describe QGP 
evolution so well. 

We thank R. Janik and P. Witaszczyk for their collaboration on related issues, A. Cher- 
nran for very valuable discussions and A. Buchel, S. Green, K. Jensen, L. Lehner, M. Lublin- 
sky, M. Martinez, W. van der Schee for helpful feedback and correspondence. This work was 
supported by the Polish National Science Centre grant 2012/07/B/ST2/03794. Research at 
Perimeter Institute is supported by the Government of Canada through Industry Canada 
and by the Province of Ontario through the Ministry of Research & Innovation. 
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Supplemental material: cancellation of ambiguities 

As stated in the section Resurgence of the Letter, the key observation is that the ambiguity 
at the leading order of the transseries is proportional to 0, so it can only be canceled by 
terms of order m — 1 or higher. The contributions to the transseries at orders 0 — 2 can be 
parameterized as 

if = l + S» ( 1 0) + + ..., 

Ir ] — c 2 ^ 2 ( a 2 ,o + •••)) (23) 

where bn ’'’ are constants, a 2 ,o = — | + and all the terms polynomial in w have been 
dropped. The constants capture the leading contributions of the ambiguous terms. 
Their values can be obtained as described in the section Hydrodynamic Gradient Expansion 
at High Orders in the Letter. The result is 

b [ 0) = e i ^r('y+l)hf ) {Co)i 

b (0) = e 2 i7r7Pr(27 + l)/ l (0)(2^ 0 ) j 

6 ( 1 ) = e -7P r(7 + 1)h W (£ o) . (24) 

The odd integer p appearing in the phase factor could a priori be different in each function, 
but for a cancellation of the ambiguity to be possible they need to be equal as written above. 
Note also that the superscript indicates the order in the transseries, so for example /i[ 0 '' is 
what in the section Hydrodynamic Gradient Expansion at High Orders section was denoted 
by hk- As stated in the Letter, the cancellation fixes the constant c up to a real number shift. 
Alternatively, one may say that the cancellation of ambiguities determines the imaginary 
part of c. 

Once the imaginary part of c is determined, conditions at higher orders in O should be 
satisfied automatically. At second order one finds 


2a 2 ,o/4° ) (£o)-M 1) (£o)=0, 

T(2 7 + l)/4 0) (2£ 0 ) - a 2 i 0 r ( 7 + l) 2 /ri 0) (£o) 2 = 0. 


( 25 ) 
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Checking this explicitly is not easy, as it requires a very accurate numerical calculation 
of the numbers h^\c, 0 ), h[ L} (C, 0 ) and hf\2£ 0 ). We founcf] that the first of the consistency 


conditions Eq. (25) is satisfied at the level of 5%. The second condition in Eq. (25) is 
much harder to check given the level of accuracy attainable using Pade approximants in the 
vicinity of branch cut singularities. 

Note finally that the term involving htp would be absent if we did not include the singu¬ 
larity at £ = 2£o in Eq. (19) in the Letter, but its presence is necessary for the cancellation 
mechanism outlined here. We expect a similar mechanism to operate at higher orders. These 
expectations are also supported by the success of the resummation described in section Re¬ 
summation in the Letter. It is also reassuring that the value of c determined there matches 
the value given in Eq. (23) in the Letter. 


* 

t 

[1] E. Shuryak, “Heavy Ion Collisions: Achievements and Challenges,” arXiv: 1412.8393 
[hep-ph]. 

[2] T. Damour, “Black Hole Eddy Currents,” Phys.Rev. D18 (1978) 3598-3604, 

[3] S. Bhattacharyya, V. E. Hubeny, S. Minwalla, and M. Rangamani, “Nonlinear Fluid 
Dynamics from Gravity,” JHEP 0802 (2008) 045, arXiv:0712.2456 [hep-th] , 

[4] N. Andersson and G. Comer, “Relativistic fluid dynamics: Physics for many different 
scales,” Living Rev.Rel. 10 (2007) 1, arXiv: gr-qc/0605010 [gr-qc] , 

[5] R. Baier, P. Romatschke, D. T. Son, A. O. Starinets, and M. A. Stephanov, “Relativistic 
viscous hydrodynamics, conformal invariance, and holography,” JHEP 0804 (2008) 100, 
arXiv:0712.2451 [hep-th]. 

[6] M. P. Heller, R. A. Janik, and P. Witaszczyk, “Hydrodynamic Gradient Expansion in Gauge 
Theory Plasmas,” Phys.Rev.Lett. 110 (2013) no. 21, 211602, arXiv: 1302.0697 [hep-th] , 

[7] I. Muller, “Zum Paradoxon der Warmeleitungstheorie,” Z.Phys. 198 (1967) 329-344. 


mheller@perimeterinstitute.ca 
michal.spalinski@ncbj .gov.pl 


1 For technical reasons, we set C\ 1 = 0 while performing the check. 









15 


[8] W. Israel and J. Stewart, “Transient relativistic thermodynamics and kinetic theory,” 
Annals Phys. 118 (1979) 341-372. 

[9] G. V. Dunne, “Resurgence and Trans-series in Quantum Theories,” Lectures given at the 
Schladming Winter School 2015 . 

[10] L. D. Landau and E. M. Lifshitz, Fluid Mechanics, Second Edition: Volume 6 (Course of 
Theoretical Physics). 1987. 

[11] W. A. Hiscock and L. Lindblom, “Generic instabilities in first-order dissipative relativistic 
fluid theories,” Phys.Rev. D31 (1985) 725-733, 

[12] P. Kostadt and M. Liu, “Causality and stability of the relativistic diffusion equation,” 

Phys.Rev. D62 (2000) 023003, arXiv: cond-mat/0010276, 

[13] R. P. Geroch, “Relativistic theories of dissipative fluids,” J.Math.Phys. 36 (1995) 4226. 

[14] R. P. Geroch, “On hyperbolic ’theories’ of relativistic dissipative fluids,” 
arXiv:gr-qc/0103112 [gr-qc], 

[15] M. Luzum and P. Romatschke, “Conformal Relativistic Viscous Hydrodynamics: 
Applications to RHIC results at s(NN)**(l/2) = 200-GeV,” Phys.Rev. C78 (2008) 034915, 

arXiv:0804.4015 [nucl-th]. 

[16] J. Bjorken, “Highly Relativistic Nucleus-Nucleus Collisions: The Central Rapidity Region,” 
Phys.Rev. D27 (1983) 140-151, 

[17] M. P. Heller, R. A. Janik, and P. Witaszczyk, “The characteristics of thermalization of 
boost-invariant plasma from holography,” Phys.Rev.Lett. 108 (2012) 201602, 

arXiv:1103.3452 [hep-th], 

[18] W. Broniowski, M. Chojnacki, W. Florkowski, and A. Kisiel, “Uniform Description of Soft 
Observables in Heavy-Ion Collisions at s(NN)**(l/2) = 200 GeV**2,” Phys.Rev.Lett. 101 
(2008) 022301, arXiv:0801.4361 [nucl-th], 

[19] M. P. Heller, R. A. Janik, M. Spaliriski, and P. Witaszczyk, “Coupling hydrodynamics to 
nonequilibrium degrees of freedom in strongly interacting quark-gluon plasma,” 

Phys.Rev.Lett. 113 (2014) no. 26, 261601, arXiv: 1409.5087 [hep-th], 

[20] A. R. Liddle, P. Parsons, and J. D. Barrow, “Formalizing the slow roll approximation in 
inflation,” Phys.Rev. D50 (1994) 7222-7232, arXiv:astro-ph/9408015 [astro-ph] . 












16 


[21] D. L. Hunter and G. A. Baker, “Methods of series analysis, i. comparison of current methods 
used in the theory of critical phenomena,” Phys. Rev. B 7 (1973) 3346-3376. 

[22] H. S. Yamada and K. S. Ikeda, “A Numerical Test of Pade Approximation for Some 
Functions with singularity,” arXiv: 1308.4453 [math-ph]. 

[23] Z. Ambrozinski and J. Wosiek, “Resummation of not summable series,” arXiv: 1210.3554 
[quant-ph]. 

[24] D. Dorigoni, “An Introduction to Resurgence, Trans-Series and Alien Calculus,” 

arXiv:1411.3585 [hep-th]. 

[25] I. Aniceto, R. Schiappa, and M. Vonk, “The Resurgence of Instantons in String Theory,” 
Commun.Num.Theor.Phys. 6 (2012) 339-496, arXiv: 1106.5922 [hep-th], 

[26] P. Argyres and M. Unsal, “A semiclassical realization of infrared renormalons,” 

Phys.Rev.Lett. 109 (2012) 121601, arXiv: 1204.1661 [hep-th], 

[27] P. C. Argyres and M. Unsal, “The semi-classical expansion and resurgence in gauge theories: 
new perturbative, instanton, bion, and renormalon effects,” JHEP 1208 (2012) 063, 
arXiv:1206.1890 [hep-th], 

[28] G. V. Dunne and M. Unsal, “Resurgence and Trans-series in Quantum Field Theory: The 
CP(N-l) Model,” JHEP 1211 (2012) 170, arXiv: 1210.2423 [hep-th], 

[29] I. Aniceto and R. Schiappa, “Nonperturbative Ambiguities and the Reality of Resurgent 
Transseries,” Commun.Math.Phys. 335 (2015) no. 1, 183-245, arXiv: 1308.1115 [hep-th], 

[30] A. Cherman, D. Dorigoni, and M. Unsal, “Decoding perturbation theory using resurgence: 
Stokes phenomena, new saddle points and Lefschetz thimbles,” arXiv: 1403.1277 [hep-th], 

[31] M. Lublinsky and E. Shuryak, “How much entropy is produced in strongly coupled 
Quark-Gluon Plasma (sQGP) by dissipative effects?,” Phys.Rev. C76 (2007) 021901, 
arXiv:0704.1647 [hep-ph], 





